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Abstract 



Using the Nikiforov-Uvarov method, the bound state energy eigenvalues and 
eigenfunctions of the ■pT-/non-'PT-symmetric and non-Hermitian general- 
ized Woods-Saxon (WS) potential with the real and complex-valued energy 
levels are obtained in terms of the Jacobi polynomials. According to the 
"PT-symmetric quantum mechanics, we exactly solved the time-independent 
Schrodinger equation with same potential for the s-states and also for any 
Z-state as well. It is shown that the results are in good agreement with the 
ones obtained before. 
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I. INTRODUCTION 

A large variety of potentials with the real or complex form are encountered in various 
fields of physics. A consistent physical theory of quantum mechanics in terms of Hermitian 
Hamiltonians is built on a complex Hamiltonian that is not Hermitian, but the energy levels 
are real and positive as a consequence of T'T— symmetry. By definition, a T'T— symmetric 
Hamiltonian H satisfies [VT, H] = 0, i.e., VTH{VT)~^ = {VT)~^HVT = H, where V and 
T are, respectively, the operators of parity (or space) and time-reversal transformations. 
These are defined according to V xV — —x, V pV — TpT — —p and TUT — —il where 
x,p, and /, are, respectively, the position, momentum and identity operatos acting on the 
Hilbert space H — L'^(Tl) and i — -s/— T. Note that this apphes only for the system whose 
classical position x and momentum p are real [1,2]. It is also known that T'T— symmetry 
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does not necessarily lead to completely real spectrum, and there are several potentials where 
part or all of the energy spectrums are complex. The Schrodinger equation (SE) for the 
real (Hermitian) potentials are investigated to generate the real energy eigenvalues which 
are of much interest [1]. Bender et al. for the first time and latter others have investigated 
several complex potentials on the T'T— symmetric quantum mechanics. The main reason for 
the growing recent interest in PT— symmetry. [1] is that the eigenvalues (spectrum) of every 
Hamiltonian is real (T'T— symmetry is exact) or come in complex conjugate pairs of complex 
eigenvalues (T'T— symmetry is spontaneously broken) [1,2,3,4,5,6,7,8,9]. Afterwards, non- 
Hermitian Hamiltonians with real or complex spectra have been studied by using numerical 
and analytical techniques [10,11,12,13,14]. 

Various different techniques have been employed in solving the above mentioned potential 
cases. One of these methods which makes it possible to present the theory of special functions 
by starting from a differential equation has been developed by Nikiforov and Uvarov (NU) 
method [15]. This method is based on solving the time-independent SE by reducing it into 
a generalized equation of hypergeometric form. Exact solution of SE for central potentials 
has generated much interest in recent years. So far, these potentials are the parabolic type 
potential [16], the Eckart potential [17,18], the Fermi-step potential [17,18], the Rosen-Morse 
potential [19], the Ginocchio barrier [20], the Scarf barriers [21], the Morse potential [22] and 
a potential which interpolates between Morse and Eckart barriers [23] . Many authors have 
studied on exponential type potentials [12,24,25,26,27] and quasi exactly solvable quadratic 
potentials [28,29,30]. In addition, Dirac, Klein-Gordon, and Duffin-Kemmer-Petiau equa- 
tions for a Coulomb type potential are solved by using this method [31,32]. The exact 
solutions for these models have been obtained analytically. 

So far, we have solved the nonrelativistic and semi-relativistic wave equations using the 
statistical model [33] , a different approach to the shifted large 1 /N expansion technique [34] 
and also the shifted large l/I expansion technique [35] with a wide group of static phe- 
nomenological and QCD-motivated potentials to produce the heavy and light quarkonium 
spectra. In addition, the energy eigenvalues of the bound states of an electron in the gen- 
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eral exponential cosine screened Coulomb potential are obtained using the shifted large 1 /N 
expansion method [36]. In this work, we solve the SE using NU method with some well- 
known WS potentials [4,13,14,37]. This potential form has been used widely in analysis of 
heavy- ion reaction and has enjoyed success [37]. It is selected for a shell model which can 
be used for describing metallic clusters in a successful way and for lighting the central part 
of the interaction neutron with one heavy nucleus [38,39]. 

In the present work, the energy eigenvalues and eigenfunctions of the Hermitian and 
non-Hermitian form of the generalized WS potential are calculated by employing the NU 
method. 

The contents of this paper is as follows: In Section II we breifly present the general- 
ized form of Woods-Saxon potential inspired from the SUSYQM. In Section III, we present 
Nikiforov-Uvarov method and also the solution of the Schrodinger equation with Hermitian 
form of the generahzed WS potential for the I — and I ^ cases. In Section IV, the 
T'T— symmetric and non-T'T— symmetric non-Hermitian potential forms are also investi- 
gated. Shapes of the potential and energy are studied in Section V. Finally, in Section VI 
we give results and conclusions. 

II. GENERALIZED WOODS-SAXON POTENTIAL 

The motion of the free electrons which have conclusive influence on the abundance of 
metallic clusters is a vital problem in the nuclear physics. These electrons are moving in well 
defined orbitals, around the central nucleus and in a mean field potential which is produced 
by the positively charged ions and the rest of electrons. In the mean field potential, the 
details of the potential are described by free parameters such as depth, width and slope 
of the potential, which have to be fitted to experimental observations. The interactions 
between nuclei are commonly described by using a potential that consist of the Coulomb 
and the nuclear potentials. These potentials are usually taken to be of the ordinary WS form 
[40] and has been widely used in analysis of heavy-ion reactions [37]. Recently, Aral [41] has 
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introduced the deformed hyperbolic potentials using the definitions of the deformed hyper 
functions. As an example, we may choose the inter-nuclear generalized WS of spherically 
symmetric form as [42] 

yir) = -V^ 1>l (1) 



where r stands for the center-of-mass distance between the projectile nucleus and the target 
nucleus. The other parameters in the potential Rq — tqA^/^ is the radius of the corresponding 
spherical nucleus or the width of the potential, A is the target mass number, tq is the radius 
parameter, Vi is the potential depth of the Coulombic part (i.e., the first term on the right 
hand side of Eq.(l)), a is the surface diffuseness parameter which is usually adjusted to 
the experimental values of ionization energies [37] and finally V2 is an introduced parameter 
for the second part of Eq.(l) (transforms like potential barrier) [42]. Further, q is the 
deformation parameter for the strength of the exponential part other than unity which is 
arbitrarily taken to be a real constant within the potential. Its worthwhile to note that the 
spatial coordinates in the potential are not deformed and thus the potential still remains 
spherical. Therefore, we start with the original SE: 

[^ + 'yi^))^i^)-E^l:{v), (2) 

where the classical phase space is assumed to be real, i.e., x and p are the standard Hermitian 
operators representing the position and momentum of a particle of mass m. Employing the 
separation of variables 

i;{v)^-R{r)Y{e,(t>), (3) 
r 

leads to the simple radial SE, for all angular momentum states, of the type 

Now, the aim is to solve the last equation for the given potential in Eq.(l) for the / = 
and / 7^ cases using the Nikiforov-Uvarov method which will be introduced breifly and 
employed in section III. 



III. APPLYING NIKIFOROV-UVAROV METHOD TO THE SCHRODINGER 

EQUATION 

The Nikiforov-Uvarov (NU) method provides us an exact solution of Eq.(4) for certain 
kind of potentials among them the one given in Eq.(l) [15]. This method is based upon the 
solutions of general second order linear differential equation with special orthogonal functions 
[43]. For a given real or complex potentials, the PT— symmetric SE in one dimension is 
reduced to a generalized equation of hypergeometric type with an appropriate s — s{x) 
coordinate transformation. Thus, it takes the form: 

r(s) + ^^'(5) + = 0, (5) 

where a{s) and a{s) are polynomials, at most of second-degree, and f (s) is of a first-degree 
polynomial. To find a particular solution for Eq. (5) by separation of variables, we use the 
transformation given by 

^{s) = cj>{s)y{s). (6) 

This reduces SE, Eq.(5), into an equation of hypergeometric type: 

a{s)y"{s) + T{s)y'{s) + Xy{s)^0, (7) 

where 4>{s) is found to isatisfy the condition 4>'{s)/4>{s) = 7r(s)/cr(s). Further, y{s) is the 
hypergeometric type function whose polynomial solutions are given by Rodrigues relation 

where Bn is a normalizing constant and the weight function p{s) must satisfy the condition 

[15] 

{a{s)p{s))' = T{s)p{s). (9) 
The function n{s) and the parameter A required for this method are defined as 
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a'js) - t{s) 
7r(s) = ± 



' "'w-n^) y_»(,)+to(,). (10) 



and 

A = A; + 7r'(s). (11) 

Here, 7r(s) is a polynomial with the parameter s and the determination of k is the essential 
point in the calculation of 7r(s). Thus, for the determination of k, the discriminant under 
the square root is being set equal to zero and the resulting second-order polynomial has to 
be solved for its roots Hence, a new eigenvalue equation for the SE becomes 

An + nr'is) + =0, (n = 0, 1, 2, ...) (12) 

where 

r(s) =f(s) + 27r(s), (13) 

and it will have a negative derivative. Therefore, we start solving SE for the VT generahzed 
WS following the NU for / = and / ^ cases as follows: 



A. Solution for the I = case 



In order to calculate the energy eigenvalues and the corresponding eigenf unctions, the 
Hermitian real- valued potential form given by Eq.(l) is substituted into the one— dimensional 
T'T— symmetrical SE with the zero angular momentum states. 



R"{x)+ . 



2m 



R{x) = 0, 



(14) 



1 + ?e-°^ (1 + qe-^^y 

where some assignments of one-dimensional parameter x — r — Rq and a — 1/a are done. 

Now, rewriting Eq.(14) by employing the convenient transformation, s{x) — — e~"^, then 
the PT— symmetrical Hermitian one-dimensional SE becomes 



(fR{s) 1 dR{s) 2m 



s ds 



E- 



Vis 



l-qs (1 - qsf 



R{s) = 0, 



(15) 



and also introducing the given dimensionless parameters: 

finally it leads into the following simple hypergeometric form given by 
d^R{s) \-qs dR{s) 



(16) 



+ 



+ 



ds'^ s{l — qs) ds [s (1 - qs)] 



X 



(-eg^ + (3q--f)s^ + {2eq - (5) s - e\ R{s) = 0. 

(17) 



Hence, comparing the last equation with the generalized hypergeometric type, Eq.(5), we 
obtain the associated polynomials as 

f{s) = l-qs, a{s) = s{l - qs), o{s) = [-eq^ + (3q - 'i) + {2eq - (3) s - t. (18) 

When these polynomials are substituted into Eq.(lO), with a'{s) = 1 — 2^5, we obtain 

n{s) = ± -^{q^ + 4eg2 - A(5q + 47 - Akq) s"^ + A{I3 - 2tq + k) s + 4e. (19) 

Further, the discriminant of the upper expression under the square root has to be set equal 
to zero. Therefore, it becomes 



(20) 



A = [4 (/5 - 2eq + k)f - 4 x 4e (g^ + Aeq^ - A(5q + 47 - Akq) = 0. 



Solving Eq.(20) for the constant k, we get the double roots as A;_|_^_ = —/3 ± g^/e (1 + ^j. 
Thus, substituting these values for each k into Eq.(19), we obtain 



^ _g£ ^ 1 [(2xA- ^TTf ) qs - 2^] ; for fc+ = -/? + q^e (l + ^), ^^^^ 
Hence, making the following choice for the polynomial 7r(s) as 



qs 1 

TTIS) = 

^ ^ 2 2 



(22) 



gives the function: 



t(s) = 1 - 2gs - 



2V^+Jl + ^) qs-2 



(23) 



which has a negative derivative of the form r(s) — — (2 + 2-y/e + ^1 + ^) q. Thus, from 
Eq.(ll) and Eq.(12), we find 



A=-,3-|(l + 2v/J)(l + ^/r7|l). 



(24) 



and 



A„ = ^2 + 2^6 + + nq + n{n-l)q. 



(25) 



After setting A„ = A and solving for e, we find: 

1 + 2n 



en{a,q) 



l + 2n+ Jl 



47 

^.2 



(26) 



Therefore, substituting the values of e, and 7 into Eq.(26) together with the transformation 
a — 1/a, one can immediately determine the exact energy eigenvalues E„{a, q) as 



E„(a,g) ^ 



2ma? 



1 + 2n 



{n{n + 1) 



2 1 + 2n + ^1 + ^ 



2 



, < n < 00, g > 1. (27) 



Following Ref.[36], in atomic units (/i = m = c = e = l), Eq.(27) turns out to be 



E„(a,g) = 



2a2 



l + 2n (n(n + l)-^ 



2 1 + 2n + 



<?2 



, < n < 00, ? > 1. 



(28) 



The above equation indicates that we deal with a family of the generalized WS spherical 
potential. Of course it is clear that by imposing appropriate changes in the parameters a 
and Vi, the index n describes the quantization for the bound energy states. In addition, if 
the parameter V2 in Eq.(28) is adjusted to zero (for s— state), solution reduces to the form 
obtained for the standard WS potential without regarding of the centrifugal barrier potential 
and the deformation parameter q in Eq.(4).^ 



^The exact bound energy eigenvalues for the s— states, in atomic units, is E„(a = l,g — 1) 
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, which is in agreement with previous works on SUSYQM (cf. formula (40) 



in Ref.[42]). 



Let us now find the corresponding wavef unctions. Applying the NU method, the polyno- 
mial solutions of the hypergeometric function y{s) depends on the determination of weight 
function p(s) which is found to be 



Pis) = (1 - qsY'-'s'^^; ^ = 1 + ^1 + i|. (29) 

Substituting into the Rodrigues relation given in Eq.(8), the eigenfunctions are obtained in 
the following form 

where C„ stands for the normalization constant and its value is l/n\. In the limit 5 — > 1, the 
polynomial solutions of Unis) are expressed in terms of Jacobi Polynomials, which is one 
of the classical orthogonal polynomials, with weight function (1 — qs^'^s'^^ in the closed 
interval [0, 1] , yielding yn,iis) — Pl^^''^~^\l — 2s) [43]. Finally, the other part of the wave 
function in Eq.(6) (with g — > 1) is found to be 

(/)(s) = (1 - s)^s^', /X = rj/2. (31) 

Combining the Jacobi polynomials and 0(s) in Eq.(31), the s— wave functions (/ — 0) could 
be determined as 

Rnis) = Dr,S^%l - s)W2p^2yi,,-l)(i _ 3^)^ (32) 

with s — — e~"^ and D„ is a new normalization constant. 

B. Solution for the / / case 
The Hamiltonian for the generalized WS potential for the / ^ case is 

To evaluate the spectra of energy eigenvalues and eigenfunctions we introduce a new effective 
potential whose form is given as: 
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eff 



( r-Ro \ 
3 V a ^ 



l{l + l)h ^a'^ 
2m 



+ 



2 ' 



(34) 



where the second term of the last expression looks like the potential barrier term of Eq.(33). 
Thus, comparing Eq.(34) with its counterpart Eq.(l), one can arrive into the transformation 



of 1^2 



2m 



Moreover, Eq.(34) can also be rewriten as 



1 



eff 



+ 



Z(Z + l)ra 



2^,2 



2 ' 



(35) 



where a — (g + l)/i?o and a — 1/a. The lowest energy levels of the potential in Eq.(34) are 
now given 



2ma'^ 



1 + 2n {ri{n + 1) 



l + 2n + + — 



, < n < oo, q>\. (36) 



For scattering processes, it has been well accepted that the surface diffuseness parameter a 
is around 0.63 /m [37,44,45]. Much larger diffuseness parameter, ranging between 0.8 and 
1.4 jm is needed in order to fit the data [37,46]. 

IV. NON-HERMITIAN POTENTIAL FORMS 
A. T^T-symmetric and non-Hermitian generalized WS case 



For T'T— symmetric and non-Hermitian potential case, we take the potential parameters, 
in Eq.(l), V2, and q as real and a — > iaj is purely imaginary. In this case, the potential 
turns out to be 



, ( q + cos ajx — i sin aix\ f q + cosajx — ismajx 
V{x) = -Vi \ — — \ +V2' 



(37) 



1 + q"^ + 2q cos ajx J \ 1 + q"^ + 2q cos ajx J 
By substituting this potential into Eq.(14) and making similar operations in obtaining 
Eq.(28), one can easily get the energy eigenvalues, in atomic units, as 



E„{ai,q) 



a/(l + 2n) (n(n + l)aj+g^ 



l + 2n+Jl-^ 



(38) 
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and 



A 



8mV2 



8mV2 



nq + n{n — l)q. 



(39) 



(40) 



Here ai is an arbitrary real parameter and i — yf—l. Thus, by choosing the parameter a 
as purely imaginary, we find the energy eigenvalues obtained for PT— symmetric and non- 
hermitian generalized WS potential are not similar to Eq.(28). A positive energy spectra is 



obtained if and only if n< - - ^1 + ^1 - , since the energy eigen- 

values of generalized WS potential are negative.^ It appears that the eigenvalues are always 



2 {l+n)aiq 



, but can be complex for 



positive real for V2 = 0, that is, En{ai, q) — j 

2 2 

V2 > Since these one-dimensional non-Hermitian Hamiltonian were invariant under 

PT— transformation, they possessed real spectra. Thus, their real spectral properties were 
linked with their PT— symmetry. 

On the other hand, following the procedures in Section III, we obtain an eigenfunction 
for the non-Hermitian potential, Eq.(37), as: 



with E„ 



n(n+l)+24 

l+2w \ 
2 



l+2n+ 



and s 



-lajx 



(41) 



. Its clear that the eigenfunctions are 



always complex for all values of V2. Therefore, if V2 = (with g — > 1), we have a simple 
wavefunction taking the form: 



Rn{x) = D„(-l)^"e-*"^^"^(l + e-'""')Pi^^-'^\l + 2e- 



(42) 



with E„ = ^ 



2 {l+n)aj 



Once we set V2 = 0, then n < J -^i^^ — 1 which agrees with Ref. [13,14]. 
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B. Non-T^T-symmetric and non-Hermitian generalized WS case 



Another form of the potential is obtained by making the selection of the potential pa- 
rameter as pure imaginary, that is, Vi is replaced by iVu, and a is replaced by iaj but V2 
remains a pure real. In this case, the Hamiltonian is non-Hermitian and non-T'T— symmetric 
having real spectra. The potential turns out to be 

2 



V{x) = -Vu 



sin ajx + i{q + cos ajx) 



+ V2 



q + cos aix — i sm aix 



(43) 



1 + + 2g cos aix ) \ 1 -|- -|- 2g cos aix ) 

By substituting this potential into Eq.(14) and making similar operations in obtaining 
Eq.(28), one can easily get the energy eigenvalues, in atomic units, as 



E„{ai,q) = i 



(l + 2n) (^f^ + n(n + l)a,) 



l + 2n + Jl - ^ 



and 



A 



z/3-|(l + z2^/i) 1^1 + ^ 



1 - 



8mV2 



(44) 



(45) 



A„ = 2 1 + iv/e + 



nq -\- n{n — l)q. 



(46) 



n^a]q^ 

The last case has real plus imaginary energy spectra. When we consider the 
real part of energy eigenvalues an acceptable result is obtained when n < 



( ilii — 1|) — — J^2^^2 1 condition. However, the energy spectrum is 



not seen at the imaginary part of energy eigenvalues, since it is independent of 



Once we set V2 = 0, then n < - 1 which agrees with Ref.[13,14]. 
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V. POTENTIAL AND ENERGY SHAPES 



It is illustrated in Figure 1, the shape of the generalized WS potential given by Eq.(l) for 
various q = 1,3 and 7 values. This is done by fixing the potential parameters Vi = 50 MeV, 
V2 = 10 MeV, ro = 1.285 fm, ^ = 56 and a = 0.65 fm [37,44,45]. The empirical values 
found by Perey et al. [47] Tq = 1.285 fm, a — 0.65 fm and A — 56 which is the geometric 
average of the target nucleus mass number 44 < A < 72 [48] are used to draw the potential 
in Figure 2 for various values of V2 = 10, 50, 100 MeV. Further, Figure 3 shows the energy 
eigenvalues in atomic units as a function of the discrete level n for various a — 0.65, 0.85, 
1.05 fm values. Of course, it is clear that by imposing appropriate changes in the parameter 
a, the index n describes the quantization of the bound states and energy spectrum. Some 
of the initial energy levels for undeformed case {q — 1) value, are presented by choosing 
Vi = 50 MeV, V2 = 10 MeV in Figure 4. Its worthwhile to note that Figure 3 and Figure 4 
have no lower bound on the spectrum for both V2 — and V2 ^ cases, respectively. In this 
regard, it is found that for V2 — case, the energy levels go into lower bound for < n < 6 
states and into higher bound for n > 6 states. However, for V2 7^ case, the energy levels 
go into lower bound for < n < 2 states and higher bound for n > 2 states. 



VI. RESULTS AND CONCLUSIONS 

In this article we used NU method and solved the radial SE for the generalized WS 
potential with the angular momentum / — (V2 — 0) and / 7^ (V2 7^ 0). A particularly 
interesting result of our investigation is that all the T'T— symmetric Hamiltonians with 
potential parameters remain all purely real have a real bound energies En with n > 
for Hermitian case and real positive in contrary to expectation if one lets a iai in 
the generalized WS potential. Therefore, for non-Hermitian case, the spectrum is real for 
V2 = but complex conjugate for some values oi V2 ^ in the generalized WS potential. 
Further, when a and Vi parameters are purely complex, it is seen that the number of 
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discrete levels for bound states is given only by the real part of energy eigenvalues. Thus, 
for a T'T— symmetric Hamiltonians the exactness of PT— symmetry implies the reality of 
spectrum. More specifically, if an eigenfunction Rn,i{s) is a PT— invariant, VT Rn,i{s) = 
Rn,i{s), then the corresponding eigenvalue of En,i is real. The exact PT— symmetry is a 
sufficient condition. But for a given T'T— symmetric Hamiltonian, it is not easy to determine 
the exactness of PT— symmetry without actually solving the corresponding radial SE. In 
this regard, the wave functions are physical and energy eigenvalues are in good agreement 
with the results obtained by the other methods [42] . 

On the other hand, the effect of the centrifugal barrier potential which goes as was 
replaced by a term having exactly the original WS form but of a second degree in order to 
reproduce its effect. This I ^ term has its physical basis arising from the superpotential 
partner of the WS potential in SUSYQM [42]. Hence, this new barrier term retakes the 
exact form of the original potential but with a small perturbed strength factor; say V2 
[42,49]. According to the complex quantum mechanics [50], the eigenvalues of the conversion 
a —>■ iaj are not simultanously eigenstates of PT— operator. 

Finally, we point out that the exact results obtained for the generalized WS potential 
may have some interesting applications in the study of different quantum mechanical systems 
and nuclear scattering. 
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FIGURES 

FIG. 1. Variation of the generalized WS potential as a function r. The curves are shown for 

various values of the deformation parameter q 



FIG. 2. Variation of the generalized WS potential as a function r. The curves are shown for 
various values of the perturbed parameter ¥2- 

FIG. 3. The variation of the energy eigenvalues with respect to the quantum number n with 
V\ = 50 MeV and V2 = 0. The curves shown are for various values of the surface diffuseness 
parameter a. 

FIG. 4. The variation of the energy eigenvalues with respect to the quantum number n with 
V\ = 50 MeV and V2 = 10 MeV. The curves shown are for various values of the surface diffuseness 
parameter a. 
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